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Abstract 

Results of a systematic study of fully integrated particle multi- 
plicities in central Au-Au and Pb-Pb collisions at beam momenta 
of 1.7A GeV, 11.6A GeV (Au-Au) and 158A GeV (Pb-Pb) using a 
statistical-thermal model are presented. The close similarity of the 
colliding systems makes it possible to study heavy ion collisions un- 
der definite initial conditions over a range of centre-of-mass energies 
covering more than one order of magnitude. We conclude that a 
thermal model description of particle multiplicities, with additional 
strangeness suppression, is possible for each energy. The degree of 
chemical equilibrium of strange particles and the relative production 
of strange quarks with respect to u and d quarks are higher than in 
e'''e~, pp and pp collisions at comparable and even at lower energies. 
The average energy per hadron in the comoving frame is always close 
to 1 GeV per hadron despite the fact that the energy varies more than 
10-fold. 



1 Introduction 

It is becoming more and more clear that results from relativistic heavy ion col- 
lisions at many different energies |l|] show striking common traits. Statistical- 
thermal models are able to reproduce particle multiplicities in a satisfactory 
manner by using a very small number of parameters: temperature, volume, 
baryon chemical potential and a possible strange-quark suppression parame- 
ter, 7s [0. We report here the results^ of an analysis of data from collisions 
at several different energies, with emphasis on the similarity of the colliding 
system. We have focussed our attention on central Au-Au collisions at beam 
momenta of 1.7A GeV (SIS) [|, ll.QA GeV (AGS) g and on central Pb-Pb 
collisions at 158/1 GeV (SPS) beam momentum [Q. As far as the choice of 
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data (and, consequently, colliding system) is concerned, our leading rule is 
the availability of full phase space integrated multiplicity measurements be- 
cause a pure statistical-thermal model analysis of particle yields, without any 
consideration of dynamical effects, may apply only in this case 0. Such data, 
however, exist only in a few cases and whenever legitimate we have extrapo- 
lated spectra measured in a limited rapidity window to full phase space. The 
use of extrapolations is more correct than using data over limited intervals of 
rapidity, especially in the framework of a purely statistical-thermal analysis 
without a dynamical model. Moreover, the usually employed requirement 
of zero strangeness {S = 0) demands fully integrated multiplicities because 
strangeness does not need to vanish in a limited region of phase space. 

In order to assess the consistency of the results obtained, we have per- 
formed the statistical-thermal model analysis by using two completely in- 
dependent numerical algorithms whose outcomes turned out to be in close 
agreement throughout. Similar analyses have been recently made by other 
authors (see e.g. [^0); however, both the model and the used data set 
differ in several important details, such as the assumption of full or partial 
equilibrium for some quark flavours, the number of included resonances, the 
treatment of resonance widths, inclusion or not of excluded volume correc- 
tions, treatment of flow, corrections due to limited rapidity windows etc. 
Because of these differences it is difficult to trace the origin of discrepan- 
cies between different results. We hope that the present analysis, covering a 
wide range of beam energies using a consistent treatment, will make it eas- 
ier to appreciate the energy dependence of the various parameters such as 
temperature and chemical potential. 



2 Data set and model description 

As emphasized in the introduction, in the present analysis we use the most 
recent available data, concentrating on fully integrated particle yields and 
discarding data that have been obtained in limited kinematic windows. We 
have derived integrated multiplicities of vr"*", A and proton in Au-Au collisions 



at AGS by extrapolating published rapidity distributions |]T2|, |13|, |I4| with 
constrained mid-rapidity value (?/nn=1-6). For proton and A we have fitted 
the data to Gaussian distributions, whilst for tt'^ we have used a symmetric 
flat distribution at midrapidity with Gaussian-shaped wings on either side; 
the point at which the Gaussian wing and the plateau connect is a free 
parameter of the fit. The fits yielded very good x^'s/dof: 0.27, 1.24 and 
1.00 for 7r+, proton and A respectively. The integrated multiplicities have 
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been taken as the area under the fitted distribution between the minimal 
Umin and maximal ?/max values of rapidities for the reactions NN — > vrNN, 
NN AK for pions and A's respectively; the difference between these areas 
and the total area has been taken as an additional systematic error. The 
area between y^^-^ and i/max amounts to practically 100% of the total area for 
pions and about 95% for A's. 

We have not included data on deuteron production because of the possible 
inclusion of fragments in the measured yields. This is particularly dangerous 
at low (SIS) energies where inclusion or not of deuterons modifies thermody- 
namic quantities like e/n [0. 

The data analysis has been performed within an ideal hadron gas grand- 
canonical framework supplemented with strange quark fugacity 7s . In this 
approach, the overall average multiplicities of hadrons and hadronic reso- 
nances are determined by an integral over a statistical distribution: 

{n,) = (2J, + 1)-^ / d'p ^— ^ — — - (1) 

(27r)^ J % ^» exp [{E, - /x ■ qi)/T\ ± 1 

where is a three-dimensional vector with electric charge, baryon num- 
ber and strangeness of hadron i as components; /j, the vector of relevant 
chemical potentials; Jj the spin of hadron i and Sj the number of valence 
strange quarks in it; the + sign in the denominator is relevant for fermions, 
the — for bosons. This formula holds in case of many different statistical- 
thermal systems (i.e. clusters or fireballs) having common temperature and 
7s but different arbitrary momenta, provided that the probability of realizing 
a given distribution of quantum numbers among them follows a statistical 



rule p, |T6[. In this case V must be understood as the sum of all cluster 
volumes measured in their own rest frame. Furthermore, since both volume 
and participant nucleons may fluctuate on an event by event basis, V and ^t 
(and maybe T) in Eq. (|1]) should be considered as average quantities ||^. 
The overall abundance of a hadron of type i to be compared with experi- 
mental data is determined by the sum of Eq. (|I|) and the contribution from 
decays of heavier hadrons and resonances: 

rii = nf™'^''^ + Br(j i)nj (2) 
j 

where the branching ratios Br(j i) have been taken from the 1998 issue 
of the Particle Data Table ||17|| . 

It must be stressed that the unstable hadrons contributing to the sum in 
Eq. (^ may differ according to the particular experimental definition. This 
is a major point in the analysis procedure because quoted experimental multi- 
plicities may or may not include contributions from weak decays of hyperons 
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and K^. We have included all weak decay products in our computed mul- 
tiplicities except in Pb-Pb collisions on the basis of relevant statements in 
ref. and about antiproton production in refs. [ITT], 113. It must be noted 



that switching this assumption in Au-Au at SIS and AGS does not affect 
significantly the resulting fit parameters. 

The overall multiplicities of hadrons depend on several unknown parameters 
(see Eq. (P) which are determined by a fit to the data. The free parameters 
in the fit are T, V, 7s and (the baryon chemical potential) whereas /ig 
and fiQ, i.e. the strangeness and electric chemical potentials, are determined 
by using the constraint of overall vanishing strangeness and forcing the ratio 
between net electric charge and net baryon number QZ-B to be equal to the 
ratio between participant protons and nucleons. The latter is assumed to be 
Z/A of the colliding nucleus in Au-Au and Pb-Pb. 

For SIS Au-Au data we have required the exact conservation of strangeness 
instead of using a strangeness chemical potential. This gives rise to slightly 
more complex calculations which are necessary owing to the very small 
strange particle production (Au-Au). The difference between these strangeness- 
canonical and pure grand-canonical calculations of multiplicities of K and A 
for the final set of thermal parameters (see Table 1) turns out to be as large 

factor 15 in Au-Au at 1.7^1 GeV. 
Owing to few available data points in SIS Au-Au collisions, we have not fitted 
the volume V nor the 7s therein. The volume has been assumed to be Airr^/S 
where r = 7 fm (approximately the radius of a Au nucleus) while 7s has been 
set to 1, the expected value for a completely equilibrated hadron gas. Since 
we have performed a strangeness-canonical calculation here, the yield ratios 
involving strange particle are not independent of the chosen volume value as 
in the grand-canonical framework. Thus, in this particular case, V is meant 
to be the volume within which strangeness is conserved (i.e. vanishing) and 
not the global volume defining overall particle multiplicities as in Eq. (p. 
Also, in order to test the dependence of this assumption on our results, we 
have repeated the fit by varying by a factor 2 and 0.5 in turn. 
A major problem in Eq. (|^) is where to stop the summation over hadronic 
states. Indeed, as mass increases, our knowledge of the hadronic spectrum 
becomes less accurate; starting from ^ 1.7 GeV many states are possibly 
missing, masses and widths are not well determined and so are the branch- 
ing ratios. For this reason, it is unavoidable that a cut-off on hadronic states 
be introduced in Eq. (^. If the calculations are sensitive to the value of this 
cut-off, then the reliability of results is questionable. We have performed all 
our calculations with two cut-offs, one at around 1.8 GeV (in the analysis 
algorithm A) and the other one at 2.4 GeV (in the analysis algorithm B). The 
contribution of missing heavy resonances is expected to be very important 
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for temperatures > 200 MeV making thermal models inherently unreliable 
above this temperature. 

Table 1: Summary of fit results. Free fit parameters are quoted along with 
resulting minimum x^'s and Ag parameters. 





Average 


Au-Au 1.7 A GeV 


T (MeV) 
^^B (MeV) 

Ts 

As 


49.6±2.5 
813±23 
1 (fixed) 
1437 (fixed) 
0.0054±0.0035 


Au-Au 11.6^ 


GeV 


T (MeV) 
liB (MeV) 

Ts 

VT3exp(-0.7GeV/T) 
As 


119.8±8.3 
553.5±16 
0.720±0.097 
2.03±0.34 
0.43±0.10 


Pb-Pb 158 A GeV 


T (MeV) 
liB (MeV) 

Ts 

V^T3exp(-0.7GeV/T) 
As 


158.1±3.2 
238±13 
0.789±0.052 

21.7±2.6 
0.447±0.025 



3 Results 

As mentioned in the introduction, we have performed two analyses (A and 
B) by using completely independent algorithms. 

In the analysis A all light-fiavoured resonances up to 1.8 GeV have been 
included. The production of neutral hadrons with a fraction / of ss content 
has been suppressed by a factor (1 — /) + /7^. In the analysis B the mass 
cut-off has been pushed to 2.4 GeV and neutral hadrons with a fraction / of 
ss content have been suppressed by a factor 7^-^. Both algorithms use masses, 
widths and branching ratios of hadrons taken from the 1998 issue of Particle 



Data Table [17 . However, it must be noted that differences between the two 



analyses exist in dealing with poorly known heavy resonance parameters. 
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such as assumed central values of mass and width, where the Particle Data 
Table itself gives only a rough estimate. Moreover, the two analyses differ 
by the treatment of mass windows within which the relativistic Breit-Wigner 
distribution is integrated. 

A summary of the final results is shown in Fig. |l|. For each analysis an esti- 
mate of systematic errors on fit parameters have been obtained by repeating 
the fit 

• assuming vanishing widths for all resonances 

• varying the mass cut-off to 1.7 in analysis A and to 1.8 in analysis B 

• for Au-Au at 1.7 A GeV, the volume V has been varied to V/2 and to 
2V (see discussion in Sect. 2) 

The differences between new fitted parameters and main parameters have 
been conservatively taken as uncorrelated systematic errors to be added in 
quadrature for each variation (see Table 1). The effect of errors on masses, 
widths and branching ratios of inserted hadrons has been studied in analysis 
A according to the procedure described in ref. |^ and found to be negligible. 
Finally, the results of the two analyses have been averaged according to a 
method suggested in ref. |5^, well suited for strongly correlated measure- 
ments. 



4 Discussion and conclusions 

From the results obtained, an indication emerges that a statistical-thermal 
description of multiplicities in a wide range of heavy ion collisions is indeed 
possible to a satisfactory degree of accuracy, for beam momenta ranging from 
1.7 A GeV to 158^4 GeV per nucleon. Furthermore, the fitted parameters 
show a remarkably smooth and consistent dependence as a function of centre- 
of-mass energy. 

The temperature varies considerably between the lowest and the highest 
beam energy, namely, between 50 MeV at SIS and 160 MeV at SPS. Sim- 
ilarly, the baryon chemical potential changes appreciably, decreasing from 
about 820 MeV at SIS to about 240 MeV at SPS. However, since the changes 
in temperature and chemical potential are opposite, the resulting energy per 
particle shows little variation and remains practically constant at about 1 
GeV per particle; this is shown in Fig. 

The supplementary 7s factor, measuring the deviation from a completely 
equilibrated hadron gas, is around 0.7 - 0.8 at all energies where it has been 
considered a free fit parameter. At the presently found level of accuracy, a 
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fully equilibrated hadron gas (i.e. 7s =1) cannot be ruled out in all exam- 
ined collisions except in Pb-Pb, where 7s deviates from 1 by more than 4a. 
This result does not agree with a recent similar analysis of Pb-Pb data |^ 
imposing a full strangeness equilibrium. The main reason of this discrepancy 
is to be found in the different data set used; whilst in ref. |§] measurements 
in different limited rapidity intervals have been collected, we have used only 
particle yields extrapolated to full phase space. The temperature values that 
we have found essentially agree with previous analyses in Au-Au collisions 
2g and estimates 11.7 A GeV 



The T value in Pb-Pb is strongly affected by high mass particle measure- 
ments, such as and S. A recent significant lowering of the S yield measured 
by NA49 |^ with respect to a previous measurement |^ results in a de- 
crease of estimated temperature value from about 180 MeV to the actual 160 
MeV. 

Forthcoming lower energy Pb-Pb and high energy Au-Au data at RHIC 
should allow to clarify the behaviour of strangeness production in heavy ion 
collision. 
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Figure 1: Comparison between particle multiplicities fitted using the thermal 
model and experimental results. 
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Figure 2: Fitted temperatures and baryon-chemical potentials plotted along 
with curves of constant energy per hadron. 
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